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Abstract. 

The origin of cosmic rays holds still many mysteries hundred years after they were 
first discovered. Supernova remnants have for long been the most likely sources of 
Galactic cosmic rays. I discuss here some recent evidence that suggests that supernova 
remnants can indeed efficiently accelerate cosmic rays. For this conference devoted 
to the Astronomical Institute Utrecht I put the emphasis on work that was done in my 
group, but placed in a broader context: efficient cosmic-ray acceleration and the im- 
plications for cosmic -ray escape, synchrotron radiation and the evidence for magnetic- 
field amplification, potential X-ray synchrotron emission from cosmic-ray precursors, 
and I conclude with the implications of cosmic -ray escape for a Type la remnant like 
Tycho and a core-collapse remnant like Cas A. 

The Astronomical Institute Utrecht had a 370 yr long history. But its successes 
during the more recent past owe much to Marcel Minnaert. In this context one some- 
times refers to the "Minnaert-school" within Dutch astronomy. Its is characterized by 
its focus on spectroscopy and knowledge of microphysical processes for understanding 
astrophysical sources. 

Here I discuss cosmic -ray acceleration by supernova remnants (SNRs), a topic that 
I started working on as PhD student at SRON, Netherlands Institute for Space Research, 
Utrecht. SRON's history is connected to Minnaert as well, as SRON was the result of 
two space science research groups in Utrecht and Leiden, both directed by two distinct 
pupils of Minnaert, respectively Kees de Jager and Henk van de Hulst. 

It is also interesting to note that the conference was held not in Utrecht, but at 
LeeuwenhorstQ Noordwijkerhout. For me this is also a special place: the high school 
I attended bears the same name and is located next to the conference center. This high 
school is remarkable for having produced two NASA Chandra fellows: Rudy Wijnands 
and myself. 



1. Introduction 

This y ear marks the 100th anniversary of detection of cosmic rays by Victor Hess (iHessI 
1912b Fl The source(s) of these highly energetic particles, and the way these particles 



1 But there is an Utrecht connection: The name derives from the abbey that was this location until the 16th 
century, which fell under the jurisdiction of the Bishop of Utrecht. 

2 See for Carlson! (120121) for the history behind the discovery of cosmic rays. 
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are accelerated, is still a matter of debate. The energy density of cosmic rays in the 
Galaxy, about 1 e Vcm~ 3 has for a long time be en attributed to the power provided 
by supernovae (e.g. iGinzburg & Syrovatskiilll964l) . Supernovae are the most energetic 
events in the Galaxy. Nevertheless, a large fraction, 10-20%, of their explosion energy 
is needed for particle acceleration, in order to explain the flux of cosmic rays observed 
on Earth. In addition, the lack of spectral features in the cosmic-ray spectrum up to 
3xl0 15 eV suggests that the sources of cosmic rays should be able to accelerate particles 
to this energyjj The power provided by supernova explosions may be used to accelerate 
particles immediately after the explosion, in the supernova remnant (SNR) phase, or 
perhaps in OB asso ciations, due to the comb ined effects of multiple supernovae and 
strong stellar winds dBvkov & Fleishmanl[l992h . 

Most evidence supports now the view that at least part of the cosmic rays originate 
from the SNR phase. Since a long time radio observations of SNRs show that relativis- 
tic electrons are present there. But over the last two decades evidence shows that SNRs 
can accelerate particles to energies of at least 10 TeV. The first proof was t he discov- 
ery of X-ray synchrotron radiation from the SN1006 ( Koyama et al.lll995h. followed 
by further evidence for hard X-ray s ynchrotron emission from Cas A dTheetalJll996t 
Allen et aT1ll997l : lFavata et al.lll997l) . Now for many young SNRs, regions close to the 



shock front have been identified whose X-ray emission is dominated by synchrotron 
radiation. Further proof for acceleration to TeV energies consists of the detection of 
TeV gamma-ray emission from many young SNRs by Cherenkov telescopes such as 
HEGRA, HESS, MAGIC and VERITAS. And since a few years the Fermi and AGILE 
y-ray observatories detected many, both young and mature, SNRs in the GeV y-ray 
rangeQ However, there is no conclusive observational evidence yet that SNRs are ca- 
pable of accelerating particles up to, or beyond, 3x 10 15 eV, or that 10% of the explosion 
energy is transferred to cosmic rays. Nevertheless, observations give us several hints 
that SNRs can indeed accelerate sufficient numbers of particles to these very high ener- 
gies. 



2. Efficient cosmic-ray acceleration 

The spectrum and flux of the cosmic rays detected on Earth require that SNRs should be 
able to transfer a relatively high fraction of their kinetic energy to cosmic rays, and that 
they are capable of accelerating to a t least 3 x 10 15 eV. A few decades this seemed im- 
possible (Lagage & Cesarsky|[l9 83) within the framework of the standard acceleration 



theory, the so-called diffusive shock acceleration (DSA) theory. According to DSA, 
particles of sufficient energy can diffusely move with respect to the overall plasma flow, 
thereby crossing the shock front. Since there is difference in plasma velocity between 
both sides of the shockjf] the particle receives a boost each time it crosses the shock. 
In fact, it can be shown that dE/E « constant. Particles advected by the shock-heated 
plasma that are too far from the shock region will be lost from the acceleration process. 



3 Galactic sources should probably even accelerate up to 3 x 10 18 eV, as only above this energy the Galaxy 
becomes transparent for cosmic rays. 

4 See lRevnoldsl Jlllh i lVmkl j2(5Tl) :lHelder et al] j20T3) for recent reviews. 

5 V2 = AV = Vi — V2 = (1 — 1/xWs, with^- the shock compression ratio and V s the shock velocity. 
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The acceleration time scale is given by approximately r acc = D/V^, with D the 
diffusion parameter given by D = 5/l m fpV part icie ~ \r\cEjeB, in which case the mean 
free path A map is assumed to be a factor 77 times the gyro-radius. For a very turbulent 
magnetic field SB/B ~ 1, 77 = 1. The length scale over which diffusive transport 
dominates over advection is given by l$g = r acc V s = D/V s . For a mean Galactic field 
of B ~ 5 /jG and V s = 5000 km , the maximum energy that can be reached in 500 yr 
is f'max < Tj~ l 6x 10 14 eV. To explain the highest energy Galactic cosmic rays, therefore, 
requires magnetic fields much higher than the mean Galactic field, and n close to one. 

If SNRs are indeed very efficient in converting kinetic energy to cosmic-ray en- 
ergy, one can no longer treat cosmic rays as test particles, but one has to consider th e 
back-reaction of the cosmic rays on the shock dynamics itself dMalkov & DrurylboOl ). 



It turns out that in general, cosmic ray acceleration can indeed be very efficient, al- 
though critical ingredient into the theory that is not well constraint is the injection of 
particles into the DSA process. Recent results show that easily more than 50% of the 



pressure can be cosmic ray pressure (e.g. iBlasi et ai1l2005l: IVladimirov et al. Il2008h . The 



result of this is that cosmic rays streaming ahead of the shock will start compression 
and heating the plasma in a so-called cosmic-ray precursor. For a normal, strong shock 
one expects a shock compression ratio of x = 4 (for y = 5/3). But the increased com- 
pression and change of velocity induced by the cosmic-ray precursors makes that for 
an efficiently accelerating shock the gas-shock compression will be x < 4. But the 
combination of shock-compression and cosmic -ray precursor compression will result 
in^ftotai > 4. The post-shock plasma temperature will be lower than expected as the gas 
is heated by a low Mach number shock. 

These effects are indicated by numerical acceleration models, but recently lVink et al 



(2010) calculated the effects just assuming thermodynamical principles alone, using a 



two component fluid: plasma and cosmic ra ys. This approach can reproduce some of 
the results of more complicated calculations. IVinketali(l2010h calculated the thermo- 
dynamics in three regions: far ahead of the shock (0), in the precursor, just before the 
plasma enters the shock (1) and in the shock-heated region (2). The results are param- 
eterized in terms of the particle pressure in cosmic-rays in region 2, w = P C r/Ptota\- It 
is assumed that the cosmic-ray pressure is constant across the shock (from 1 to 2). The 
conservation equations to be solved are just the usual shock equations for mass conser- 
vation and pressure equilibrium in each region [pv], [P + pv 2 ]0 But for the energy flux 
it assumed that some fraction, e, of the free energy can be carried away in the form of 
escaping cosmic rays: (P2 + U2 + = (Pq + uq + (1 - e)|pv^)vo. 

The physical reason that energy escape is necessary for efficient cosmic-ray accel- 
eration is that a shock implies a rapid change in velocity across a region. This can be 
either accomplished by increasing the entropy, as is the case in a strong non-cosmic-ray 
accelerating shock. Or energy has to leak out of the system, if cosmic rays are acceler- 
ated, because the cosmic-ray acceleration does not result in a big jump in entropy. 

The resulting equations for the shock compression ratio, the dependence of cosmic- 
ray pressure on cosmic-ray escape, and the post-shock temperature are graphically 
represented in Fig. [T] Noteworthy is that for infinite Mach number the relation be- 
tween partial cosmic -ray pressure on the one hand, and the overall compression ra- 
tio X\2 = Aftotai an d main shock compression ratio X2, on the other hand, is given by 
w ~ ixn -Xi)lkX\i ~ !)• Another interesting result that is usually taken from numeri- 



6 The brackets indicate conserved quantities. Velocities are evaluated in the frame of the shock. 
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Figure 1. Shock relations for a mixture of hot plasma, y — 5/3 and cosmic 
rays, y - 4/3, as a function of w, i.e. partial cosmic-ray pressure. Shown are 
the total ixn) and gas-shock (xi) compression (top, left); the energy-escape flux, 
e (top right); the attenuation of the post-shock temperature (bottom,left); and the 
ratio of the energy flux escaping and the flux advected with the plasma flow (bottom 
right). Mach numbers (M) are distinguished by diffe rent colors ( online edition only). 
Crosses indicate calculations based on the model of iBlasi et al.l d2005l) for M = 100. 
The downsloping parts should be considered unstable solutions. 



cal simulations, but can be analytical derived, is that the maximum overall compression 
ratio is reached for a compression of X2 - Tgas/(Tgas - 1) = 2.5 at the main shock. 
Finally, it turns out that for Mach numbers M < 6 efficient cosmic ray acceleration 
cannot be sustained (blue lines in Fig. [D. On face value this means that intra-cluster 
shocks cannot be efficient particle accelerators. 

The calculated plasma temperatures can be compared to measured plasma temper- 
atures to infer the cosmic -ray pressure, if the shock velocity is known. This has been 
investigated now using optical Ha spectroscopy, which indicates that indeed cosmic- 
ray pressure is contributing substantially to the overall pressure in some young SNRs, 
with w > 20% (He lder et al.l[2009Ll2010h . 



3. X-ray synchrotron radiation and evidence for magnetic field amplification 

As explained above, the requirement that SNRs should be able to accelerate above 
3x 10 15 eV in order to explain the cosmic ray spectrum on Earth, requires both magnetic 
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Figure 2. Top: X-ray image of Tycho's SNR, as observed by Chandra in 2009. 
Red indicates Fe-rich plasma, green Si-rich plasma, and blue strong continuum emis- 
sion (X-ray synchrotron radiation). Right: a detail of the Western shock region. 
Bottom panel: Emission profile from a 30° segment of the Western shock region in 
the 4.3-5.2 keV band. The blue line shows a emission from a uniformly emitting 
shell with thickness /. The red and line show models in which the emissivity falls 
of behind the shock as exp(-(>o - r)/Z), but the green line has an additional precur- 
sor component. The models take into account spherical projection and instrumental 
resolution. (See the online edition for colors.) 



fields that are larger than the mean Galactic magnetic field and that the magnetic fields 
are highly turbulent (i.e. ^* 1). 

The detection and properties of the X-ray synchrotron emission from young SNRs 
has shown that this may indeed be the case. First of all, the maximum energy that 
electrons can be accelerated to is limited by radiative (synchrotron) losses, which have 
a typical time scale of ti oss = 636 /B 2 E s. Equating T] oss with the typical accelera- 
tion time scale r acc (§12) shows that for electrons x VlKnceB). The resulting 
typical maximum photon e nergy turns out to be independent of magnetic field (e.g. 



Aharonian_&Atoyanj[l999): hv = 1AE 2 B keV « \An-\VJ50m\m^ 1 ) 2 keV. Since 



all young SNRs with V s « 2000 - 6000 km/s show evidence for X-ray emission 
(hv ~ 0.1-10 keV), the magnetic field must be rather turbulent, i.e. r] < 10. This 



solves one of the problems for cosmic -ray acceleration noted by lLagage & Cesarsky 
(119831 ). 
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ihe A-ray synchrotron emission irom many young mnks is clearly continea to re- 
gions close to the shock front (e.g. lVink & LaimngKOOSjBamba et al.l2005l : lBerezhko et al 
2003; He lder et alJl2012h . This is in particular true for Cas A, Tycho's SNR and Ke- 
pler's SNR. Fig. |2] shows the Chandra X-ray image of Tycho's SNR (SN 1572). It 
beautifully shows the hot thermal emission from supernova ejecta, shown as a fluffy 
structured colored red and green in the picture. It also shows a very narrow, purplish 
blue, filament forming the outer boundary of the SNR, which is caused by X-ray syn- 
chrotron radiation. The thinness of this region can be understood in two ways: one 
by noting that as the plasma flow is advected away from the shock with a velocity 
Av = Vs(l - l/x), carrying away cosmic -ray electrons. Because the electrons suffer 
radiative losses, they will after traveling a distance / = V2Ti oss no longer emit X-rays. 

Another way to look at it is, to say that only electrons that remain within a diffusion 
length /diff of the shock (Sect. 12) will be able to compensate radiative losses with shock 
acceleration energy gains. Combining these two ideas in fact shows that the physical 
widths of the X-ray synchrotron e mitting filaments depend on the average magnetic 



field strengths behind the shock (e.g Berezhk o & Volkll2004klVinkl l2005; P arizot et al.l 
2006l : IVinkll2012h : 



B « 26( ^_)- 2/3 ^/3(^ 4 - I)" 1/3 ^G, (1) 

V 1.0x 10 18 cmVr v ^ 

with /synchr the width of the synchrotron filaments and xa the overall compression ratio 
in units of x = 4. Interestingly, this equation is independent of V s , but it assumes that 
the synchrotron emission is observed with photon energies near the maximum energy. 

As an example, Fig. |2] (lower panel) shows the X-ray synchrotron surface bright- 
ness in the Western region. The actual width depends on how the emissivity falls off 
with distance to the shock. Acceleration theory predicts an exponential fall-off, but as 
the figure indicates this does not give a very good fit. A uniform shell with a sharp edge, 
provides a better fit, as it falls-off less sharply from the maximum to larger radii. The 
exponential model gives a width of / = 1.2" corresponding with / = 5.4 x 10 16 cm, and 
implying a magnetic field of B ~ 230 /uG. On the other hand, for a uniform emissivity, 
I = 3.3", which implies B x 100 //G. 

During the conference, Colin Norman reminded me of the effort of Bram Achter- 
berg to use radio-synchrotron radiation f rom Tycho's SNR to detect radio synchrotron 
emission from the shock precursor ( Achterberg et al.fl994 ). In principle, in X-rays this 



should be easier to do, as the shock precursor length should be larger in X-rays. On 
the other hand, the emissivity scales as S oc B a+ , with a the energy flux index. In 
X-rays a « 2, whereas in the radio a « 0.6. Since the shock compression makes that 
in the precursor the magnetic field is lower by a factor ~ 3 for a compression factor 
of 4, one does not expect to detect the X-ray precursor. But if the shock is efficiently 
accelerating, the main shock may have a compression factor of 2.5 (§0), which makes 
that a turbulent upstream magnetic field is on average only a factor 1.9 weaker, and 
the synchrotron emissivity only 14% of the post-shock emissivity. Indeed, allowing 
for precursor emission that is 14% weaker does improve the fit, and results in an even 
smaller width in the post-shock region / = 0.4" (implying B ~ 480 /iG), and a precursor 
length of 1.8". I cannot claim that this is inde ed a detection of t he precursor, but it is a 
hint that needs to be further explored (see also Bamb a et al.ll2005h . 

It turns out that applying Eq. [TJto X-ray synchrotron radiation from young SNRs 
implies that all these SNRs have magnetic fields close to the shock front that are 
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substantially large r than the mean Galactic field (iBamba et al.ll2005l : IVolk et all 12005 



Held er et al.ll2012 ). This implies that protons and other atomic nuclei can be accelerated 
to energies close to 3 x 10 15 eV. Moreover, there seems to be a trend that SNRs in lower 
density environments have lower magnetic fields. The correlation with shock velocity 
indica tes that B 2 « p vf , with B = 2-3 (IVolk et al.ll2005l : I Vinkl 120081 : iHelder et all 



20121) . This is in agreement with r ecent theor i es about magnetic field amplification due 



to cosmic ray induced turbulence (Bell 2 0041) . IHelder & Vi nk (2008) reported evidence 
for SNR Cas A that also the reverse shock, which heats the ejecta, is a dominant source 
of X-ray synchrotron radiation. Also in the reverse shock region the magnetic field 
must be 100 - 500 juG judging from filament widths. What is remarkable about this is 
that the ejecta are thought to have a low magnetic field due to the tremendous expansion 
of the material. That the magnetic field is nevertheless high suggests that only small 
seed magnetic fields need to be present to amplify the magnetic fields by large factors. 
Or perhaps one should even call it magnetic field creation. 

The ideal environment to accelerate to high energies are SNRs that evolve into 
the winds of their progenitors. In these winds the density drops of a s p oc 1 /r 2 (e.g. 



Schure et aLKOlOh . And as a result the density is high during the first century of the life 



of a SNR, when the shock velocity is also high. The flux of particles entering the shock 
that can potentially be accelerated is also large early on, since F cam oc p4nr 2 V s . This is 
^cam x 47r(M/v w )V 4 for SNRs in a dense wind. So more particles are accelerated early 
on, if SNRs evolve in a dense wind. The magnetic field amplification ensures a high 
magnetic field, and particles can be, on average, accelerated for a longer time. 



4. Cosmic-ray escape 

The idea that core collapse SNRs evolving in the dense winds of the progenitor stars 
are better at accelerating cosmic rays than Type la SNRs is strangely enough not sup- 
ported by the most recent y-ray observations. Cherenkov telescopes show that most 
young SNRs emit TeV radiation, but for many it is uncertain whether this is caused by 
electrons (through inverse Compton scattering) or cosmic-ray nuclei (pion decay). The 
GeV telescopes Fermi and AGILE have now greatly increased the number of SNR y- 
ray sources, with many of them being older, core collapse SNR s, often associated w ith 
dense environments, and spectra that cut-off above ~ 10 GeV (Helder et al.ll2012l for 



a review). Some of these sources have associated TeV sources that seemed displaced 
from the SNR itself, sugge sting that a nearby cl oud is hit by cosmic rays that have es- 



caped from the SNR (e.g. lAharonian et aUl2008h . So clearly there is strong connection 



between core-collapse supernovae and cosmic ray production. The mature SNRs seem 
to have lost the highest energy particles. 

When it comes to young SNRs, y-ray observations do not so clearly indicate that 
a sizable fraction of their energy is in the form of cosmic rays. Naively one would 
expect that a radio-bright and dense core collapse SNR like Cas A would have a large 
cosmic ray content, but Ferm i observations in dicate that less than 4% of the explosion 
energy went to cosmic rays (lAbdo et al.ll2010h . In contrast, the Type la SNR Tycho 
seems to have a put a larger fraction of its kinetic energy into accelerating cosmic rays 
(lAcciari et al.ll20li1: iGiordano et al. 120121) . Given that in § [3] I have argued that core 



collapse SNRs evolving in dense winds are probably better accelerators, this contrast 
between these two young SNRs seems odd. 
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However, one should also consider that a core collapse SNR like Cas A may have 
accelerated more cosmic rays in the first century after the explosion, whereas Tycho 
may now reach its peak, for the reasons given in § [3] The relative lack of y-raysfrom 
Cas A would then imply that most cosmic rays have escaped the SNR. As explained in 
§ |2]escape is a necessary element for efficiently accelerating shocks. 

There can be an additional reason that explains the difference between Tycho and 
Cas A. The emission in pion-decay depends on the local cosmic-ray number den- 
sity hcr, the density of the local medium, n p and on the emitting volume: L p i on oc 
J ncR n pdV. Due to diffusion, cosmic-rays will occupy a region that is larger than the 
SNR itself, say with radius Rdis- The cosmic-ray number density can be expressed 
in terms of the ratio of the total number of cosmic rays and the emitting volume 
n C r = N CT /V. For a Type la, like Tycho, which probably evolves in a medium with 
constant density, one sees that the total luminosity of the region containing Tycho only 
depends on N CI and n p , and is independent of Rug: L p i on oc J ncsJi v dV « N cr n v . This 
is in contrast to a SNR evolving in a wind: 

rR&m y[ N M 

Jo ^nr v w K diff v w 

The contrast between Tycho and Cas A could, therefore, imply that part of the 
y-ray emission of Tycho may come from a region outside the SNR, whereas for Cas A 
most y-ray emission comes from the SNR shell. Of course, if in the case of Type la 
SNRs /?diff becomes too large, the emission may fall below the surface brightness limit 
of the telescope, since the surface brightness scales with L p i on /R^ iS . In the future, with 
a telescope like CTA one may want to search for these cosmic-ray haloes around SNRs. 

With this view toward the future I like to conclude my chapter on recent research 
done by me and my research group at the Astronomical Institute Utrecht. For me this 
closes off a chapter in my research life of about six years, which amounts to about 1.6% 
of the total history of the astronomical research in Utrecht. 
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